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CONVERGENCE  PERFORMANCE  OF  ADAPTIVE  DETECTORS, 

PART  2 


1.0  INTRODUCTION 

The  general  problem  of  detecting  signals  in  a  background  of  Gaussian  noise  for  an  adaptive 
array  was  addressed  by  using  the  techniques  of  statistical  hypothesis  testing  (1).  In  Ref.  1,  the  prob¬ 
lem  is  formulated  as  a  binary  hypothesis  test  where  one  hypothesis  is  noise  only  and  the  other  is  sig¬ 
nal  plus  noise.  A  given  input  data  vector  (called  the  primary  data  vector)  is  tested  for  signal  pres¬ 
ence.  Another  set  of  signal-free  data  vectors  (called  the  secondary  data  vectors)  is  available  that 
shares  the  unknown  covariance  matrix  M  of  the  noise  in  the  primary  data  vector.  A  likelihood  ratio 
decision  rule  is  derived,  and  its  performance  is  evaluated  for  the  two  hypothesis. 

Kelly’s  detector  [1]  uses  the  maximum  likelihood  (ML)  estimates  for  the  unknown  parameters  of 
the  likelihood  ratio  test  (LRT).  The  unknown  parameters  are  the  spatial  covariance  matrix  and  the 
unknown  signal’s  complex  amplitude  (assumed  in  Kelly’s  analysis  to  be  a  nonrandom  constant).  This 
detection  scheme  is  commonly  referred  to  as  the  generalized  likelihood  ratio  test  (GLRT)  and  is 
referred  to  here  as  such. 

Convergence  results  for  PD  and  PF  are  presented  in  Ref.  1 .  Expressions  for  PD  and  Ph  are 
derived  that  are  a  function  of  the  number  of  statistically  independent  secondary  data  vectors;  the 
number  of  input  channels  (V;  the  detector  threshold;  and  the  input  signal-to-noise  (S/N)  power  ratio. 
Note  that  PF  did  not  depend  on  M  (a  statistical  measure  of  the  external  noise  environment).  Hence, 
this  detector  exhibited  the  desirable  constant-false-alarm-rate  (CFAR)  property  of  having  the  PF  be 
independent  of  the  covariance  matrix.  Additional  research  in  the  area  is  contained  in  Refs.  2  through 
5. 


A  mean  level  adaptive  detector  (MLAD)  is  a  more  easily  implementable  adaptive  detection 
scheme.  The  MLAD  is  essentially  an  adaptive  matched  filter  (AMF)  followed  by  a  mean  level  detec¬ 
tor  (MLD)  [6,7],  Input  samples  used  in  determining  the  MLD  threshold  are  derived  from  a  block  of 
data  passing  through  the  AMF.  The  squared  magnitude  of  each  of  these  same  samples  as  processed 
through  the  AMF  is  used  as  a  test  statistic  and  is  compared  against  an  MLD  threshold  (an  average  of 
the  instantaneous  powers)  that  does  not  contain  the  given  test  statistic  sample.  We  further  clarify  the 
implementation  terminology  by  calling  this  an  MLAD  with  concurrent  data  samples.  In  Part  1  of  this 
report  18),  an  analysis  is  performed  for  an  MLAD  with  nonconcurrent  data;  i.e.,  the  MLD  test 
statistic  is  derived  from  a  set  of  data  that  are  statistically  independent  of  the  test  statistic  of  the  pri¬ 
mary  data  vector. 

Here  we  compare  the  detection  performance  of  the  two  detection  schemes;  GLRT  vs  MLAD. 
For  this  analysis,  we  assume  that  the  complex,  desired  signal  amplitude  is  a  complex,  zero-mean 
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Gaussian  random  variable  (RV)  of  unknown  variance  with  independent  and  identically  distributed  real 
and  imaginary  parts  (the  magnitude  of  this  amplitude  is  Rayleigh  distributed).  Under  the  GLRT,  we 
would  have  had  to  reformulate  Kelly’s  detector  with  the  variance  of  the  unknown  signal  amplitude  as 
an  unknown  parameter  and  find  the  ML  estimate  of  this  quantity.  This  is  mathematically  tedious.  In 
lieu  of  implementing  this  new  GLRT.  we  chose  to  evaluate  Kelly’s  GLRT  as  it  is  defined  in  his 
paper.  As  noted  by  Kelly,  no  optimality  properties  are  claimed  for  this  test.  The  form  of  the  test  is, 
however,  reasonable. 

2.0  GENERALIZED  LIKELIHOOD  RATIO  TEST  (GLRT) 

2.1  Detector  Form 

Kelly  [1]  gives  a  mathematical  formulation  of  the  adaptive  detection  problem  that  leads  to  the 
GLRT.  We  now  summarize  that  formulation.  Two  sets  of  input  data— primary  and  secondary— are 
used.  We  assume  that  the  secondary  inputs  do  not  contain  the  desired  signal.  Set 

X  —  N  x  K  matrix  of  secondary  input  data.  The  nth  row  represents  the  K  samples  of  data  on 

the  nth  channel,  where  n  =  1,  2 . N.  The  samples  in  the  &th  column  are  assumed 

time-coincident,  where 

x  is  the  primary  data  vector  of  length  N  and 

s  is  the  desired  steering  vector  of  length  N. 

Consider  the  two  hypothesis 

H0  :  x  =  n 

H  |  :  x  =  n  +  a  s. 

where  H0  is  the  noise  only  hypothesis,  n  is  a  noise  vector  of  length  N,  and  H]  is  the  signal-plus- 
noise  hypothesis,  where  a  is  the  unknown,  complex  signal  amplitude.  We  assume  the  following: 

(1)  Input  noises  are  complex,  zero-mean  stationary  Gaussian  RVs.  The  real  and  imaginary 
parts  of  a  given  input  noise  sample  are  independent  and  identically  distributed  (IID)  with 
respect  to  each  other.  An  RV  with  these  characteristics  is  called  a  circular  Gaussian  pro¬ 
cess. 

(2)  Input  noise  samples  are  temporally  statistically  independent. 

(3)  The  secondary  data  is  statistically  independent  of  the  primary  data. 

(4)  The  desired  signal  is  present  in  the  primary  data  vector.  It  is  not  in  the  secondary  data. 

The  GLRT  is  formulated  as  follows.  Find  the  probability  density  function  (PDF)  under  each 
hypothesis  over  all  measured  data.  For  this  problem,  this  is  straightforward  since  the  sample  vectors 
are  assumed  to  be  independent  and  each  vector  has  an  associated  /V-dimcnsional  Gaussian  PDF.  If 
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there  are  any  unknown  parameters,  maximize  the  PDF  of  the  inputs  over  all  unknown  parameters  for 
each  of  the  two  hypothesis.  The  maximizing  parameter  values  are,  by  definition,  the  ML  estimators 
of  the  parameters.  Hence,  obtain  the  maximized  PDFs  by  replacing  the  unknown  parameters  by  their 
ML  estimates.  Find  the  ratio  of  the  resultant  maximum  of  PDFs  (the  ratio  of  the  PDF  under  H  (  to 
the  PDF  under  Ha).  Check  this  ratio  to  see  if  it  exceeds  a  preassigned  threshold  t. 

Kelly  shows  that  the  GLRT  for  the  adaptive  detection  problem  is  given  by 


SH  k;'  x| 


", 

_  > 

(sH  R~'  s)  [1  +  xH  k;'  x]  < 


/. 


where 


(1) 


Rx=XXH.  (2) 

and  H  denotes  the  conjugate  transpose  matrix  operation.  We  recognize  Rx  as  proportional  to  the  ML 
estimate  of  the  input  covariance  matrix.  We  note  also  that  the  desired  signal's  unknown  complex 
amplitude  a  has  been  estimated  and  is  accounted  for  in  Eq.  (1). 

2.2  Statistically  Equivalent  GLRT 

Here  we  derive  a  statistically  equivalent  GLRT  that  allows  us  to  formulate  simply  the  detection 
and  false  alarm  probabilities  of  the  adaptive  detector.  As  in  Refs.  1  and  9.  we  can  matrix  transform 
the  input  vectors  by  an  (V  x  N  matrix  A,  which  has  the  properties  that  the  input  noise  vectors  are 
spatially  whitened,  each  input  element  has  noise  power  normalized  to  one,  and 

/4s  =  (0,0,  •  •  •  0,  (sH  M~'  s),/2)  =  s0,  (3) 

where  all  of  the  desired  signal  has  been  placed  into  A/th  channel  (note  in  Ref.  1.  the  signal  was  placed 
into  the  first  channel;  for  analysis  purposes,  we  place  the  signal  into  the  Mh  channel). 

In  addition,  set 


z  =  Ax,  and 

(4) 

T 

II 

N 

(3) 

The  elements  of  vector  z  (under  H())  and  the  elements  of  the  vectors  representing  the  columns  of  Z 
teach  column  represents  the  transformed  secondary  data  across  the  array  at  a  given  instant  of  time), 
are  now  spatially  independent  with  each  element  having  power  equal  to  one.  As  shown  by  Ref.  I, 
the  transformed  GLRT  is  given  by 
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where 


So  /?.  1  z 


1  I 


(So  R;  S0)(l 


ZH  R.  '  2) 


W, 

> 

< 

H  n 


R.  =  ZZH. 


(6) 


(7) 


We  note  that  if  Q  is  any  K  x  K  unitary  matrix  (i.e.,  QQH  =  /*  where  IK  is  the  K  x  K  iden¬ 
tity  matrix),  then 


R.  =  ZZH  =  ZC?£?H  Zh. 


(8) 


In  the  forthcoming  development,  we  will  transform  Z  by  a  series  of  unitary  matrix  transforms  that 
zero  many  of  the  elements  of  Z.  This  development  is  similar  to  that  of  Rader  and  Steinhardt  [  10], 


Set  Z  =  (z^)  =  Z0.  We  can  show  by  construction  that  there  exists  a  K  x  K  unitary  matrix 
*o  such  that 


Vn 

0 

o 

•  •  • 

ziV 

zs  y 

<-n\ 

ztfi 

- 1 

where  is  a  function  of  only  z\k,  k  =  1,2 . K.  and 


(9) 


Vi  | 


(10) 


Note  that  after  this  transformation.  R.  is  dependent  only  on  the  RVs,  yM,  n  =  2,  . . .  ,  N, 
k  =  \ ,  2,  ....  K.  Hence,  we  have  reduced  the  total  number  of  RVs  on  which  R.  is  dependent  by 
K  -  1.  In  addition,  yn  is  statistically  independent  of  the  z(n[',  and  the  elements  of  the  second 
through  Mh  rows  are  statistically  independent  with  each  element  having  a  power  equal  to  one.  (This 
follows  because  a  vector  of  IID-complex,  circular  Gaussian  RVs  is  the  same  statistically  after 
transformation  by  a  unitary  matrix).  We  note  that  vn  has  a  x  PDF  of  the  order  2 (K  -  1)  with 
<r  -  0.5. 
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Furthermore,  there  exists  a  K  x  K  unitary  matrix  such  that 

Z2  =  Z,*,. 


(ID 


and 


1  0 
o 


yn  0  0  ...  0 

zV,’  >22  0  ...  0 

4V  Z®  ...  zS! 


(12) 


L')vi  z]d\  z'fji,  . . .  z(sk J 

where  'I'l  is  a  (AT-1)  x  (A'-l)  unitary  matrix  whose  elements  depend  only  on  C2/0 
k  =  2.  3,  ...  K,  and 


Vv>  = 


L  \&'\2 


k  =2 


1/2 


(14) 


Again,  note  that  we  have  reduced  the  number  of  RVs  on  which  R-  depends.  Also,  all  of  the  elements 
of  Z2  are  independent.  The  i  =  1,  2  are  IID  circular  Gaussian  RVs  with  power  equal  to  one 
and  v 22  is  a  x  PDF  of  the  order  2{K  -  2)  with  a2  =  0.5. 

This  procedure  can  be  reiterated  so  that  the  final,  transformed  matrix  is  given  by 


where  ZN  is  an  N  x  N  matrix  defined  as 


Z,v  = 


Z,v 

=  l^v  1  0], 

V|| 

0  0  .  .  . 

0 

4V 

v’22  0  •  •  • 

0 

-11! 

<-31 

<$2  3*33  •  •  • 

0 

-(h 
'■hi  1 

zfo  ... 

>’,v,v 

(10) 


(16) 
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We  note  that  the  N  +  1  through  Kth  columns  of  ZjV  are  zero-filled  and  are  denoted  by  0.  Hence,  R. 
only  depends  on  lheN(N  +  l)/2  RVs.  All  of  the  elements  of  ZN  are  statistically  independent.  The 

z\ll ,  /  =  1,  2 . N  -  1  are  IID  circular  Gaussian  RVs  with  power  equal  to  one  and 

y,m.n  -  1,2,...  ,  N  is  a  \  PDF  of  the  order  2 (K  —  n)  with  <r  =  0.5.  Thus, 


R-  —  Zv  Z,v  —  Z; 


>H 

'N 


(17) 


We  perform  one  last  matrix  transform  on  the  input  data.  We  multiply  z  by  an  N  x  N  unitary 
matrix  transform  B ,  such  that 


B’  0 

0  1  ’ 


B  z 


0 

0 


0 

^2 


(18) 


(19) 


and 


vi 


"a-  i 

L 

fn  =  l 


l^l2 


1/2 


(20) 


v-> 


•:.v . 


(21) 


where  B'  is  an  (N  -  1)  x  (/V  —  1)  unitary  matrix.  We  note  that  B s0  =  s0,  and  if  we  set  Zg  =  BZ, 
then 


s„  -  S(?  k:' 


s,7  R: 


So. 


and 


S„H z  =  SoH zs,  (22) 

zH  r:  =  7"  k:'  /.B. 
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where  RZh  =  ZB  ZB .  The  statistical  properties  of  ZB  are  identical  to  Z  and,  hence,  Zs  can  be 
reduced  as  given  by  Eqs.  (15)  and  (16).  We  point  out  that  v,  and  v2  are  independent  and  that  \-2  is  a 
circular  Gaussian  RV  with  its  power  depending  on  the  given  hypothesis  (under  Ht).  the  power  equals 
one),  and  v,  is  a  x  PDF  of  order  2 (N  -  1)  with  a 2  =  0.5. 

As  a  result  of  the  transformations  of  the  input  data,  we  see  that  Eq.  (6)  can  be  rewritten  as 


S„H  (Z\;  Z,v ,  '  ZB  I  2 


[S0H  (ZN  Z,v)“ 1  Sol  It  +  Zs  (Z,v  Z,v)”!  zBi 


H\ 

> 

< 

Ho 


(23) 


Define  m  j j  —  \'s- ll22  ~  }’nN’u 21  ~  '!v”,v -'i  • 


U  = 


0 

U  22 


\  =  (v i , Vi)T ,  and  I0  =  (0, 1)T.  It  is  straightforward  to  show  that  Eq.  (23)  is  equivalent  to 


i,T  (UUH)~]  v 


(l0r  (UUH)-[  10)  11  +  vH  (UUH)-'  v] 


> 

< 

Ho 


(24) 


In  going  from  Eq.  (23)  to  Eq.  (24).  we  have  taken  advantage  of  the  fact  that  if  C  is  a  nonsingu¬ 
lar  N  X  (V  lower  triangular  matrix  that  is  partitioned  as 

Cl|  C  |2 
^  Ci]  C22 

where  C22  is  a  n  X  n  matrix,  and 

fc"  Cir 

C_l  =  [c21  C22  ' 

where  C22  is  a  n  X  n  matrix,  then  C22  =  C22  ■ 

Now 

i£(c/(/Hr'  i„  =  4-.  (25) 

U  22 


iH 


(UUH) 


v  -  = 


«Ti 


4 

U  ii 


Mil  V 


11*2 


'  2 1  1  I 


(26) 
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and 


vH  (UUH)~X  v 


1 

T  *> 

U' ji  Wi 2 


I  W  1 1  2 


I  2  .  2  2 

\* 2  i  V  i  |  4-  U22  V 1 


(27) 


By  inserting  Eqs.  (23)  through  (27)  into  Eq.  (24),  we  show  that  the  GLRT  is  equivalent  to 


where 


1  2  >  /  2  .  2  v  t* 

W  1 1  V 2  —  V2  ^21  I  ^  (W  1 1  +  V 1)  M22*\ 

H0 


(28) 


(29) 


2.3  Probability  of  Detection 

Under  the  H ,  hypothesis,  we  assume  that  the  desired  signal's  amplitude  (or  magnitude)  is  Ray¬ 
leigh  distributed  and  that  the  signal's  phase  is  uniformly  distributed  between  (0,  2ir).  This  implies 
that  the  desired  signal  itself  is  a  complex  circular  Gaussian  RV;  therefore,  let  the  desired  signal's 
input  power-per-channel  before  any  matrix  transformation  be  equal  to  a~ .  After  the  A  matrix 
transformation  (whitening  and  placing  the  signal  into  the  Nth  channel),  the  signal  power  in  the  Mh 

channel  is  o;  --  as  sH  M  's.  Thus  under  Ht,  v>2  is  a  complex  circular  Gaussian  RV  with  power 
equal  to  o~  +  1 ,  where  the  1  represents  the  noise  power-level-per-channel  after  the  A  matrix  transfor¬ 
mation.  Also,  u 21  is  a  complex  circular  Gaussian  RV  with  power  equal  to  one.  We  can  rewrite  Eq. 
(28)  as 

«. 

I  a  1 2  <  T\  (30) 

Ho 

where 


a 


U  1 1  V’ t  V  |  U  2  ] 


0  I  1  (<7‘  4-  1)  +  vq  ) 


,2  >1/2 


,  and 


(31) 


+  I’T 


«ii  (<j;  +  1 )  +  v] 


U  A  A 


T. 


(32) 


It  is  straightforward  to  show  that  a,  when  conditioned  on  «n  and  V|.  is  a  complex  circular  Gaussian 
RV  with  power  equal  to  one.  It  is  well  known  (11)  that  the  conditional  probability  of  detection  is 
given  by 
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,  \  Mj|  +  V]  I 

P(D\uu,  uZ2'  =  e  =  exp  <  -  — - ; - -  uzzT  >  ,  (33) 

“h(«i  +  l)  +  vi  j 

where  P(D  |  •.  )  denotes  the  conditional  probability  of  detection. 

We  set 


03 

1! 

ioisj 

to 

(34) 

p  =  «ii- 

(35) 

v  =  V|,  and 

(36) 

^  V 

(37) 

P 

Then  Eq.  (33)  becomes 


P(D\r.r,)  =  exp  -  ,  *  +/  » ,7 

!  a-  +  1  +  r  \ 


(38) 


The  PDFs  of  rj.g,  and  N  are  x2  of  order  2 (K  -  N  +  1).  2{K  -  N  +  2).  and  2(iV 
with  cr  =  0.5,  and  are  given  by 


1).  respectively 


P  7,0?)  = 


1 


(K  -  A3) ! 


7  ^ 


Ar_A'  c-’  .17  >  0. 


(39) 


PJp)  = 


l 


(K  -  N  +  1)! 


7/'--v  +  1e-"  ,/t  >  0, 


and 


(40) 


3M")  = 


1  ,v  -2 


(A3  -  2)! 


”  ,v  >  0. 


(41) 


By  using  the  elementary  probability  theory,  it  is  straightforward  to  show  that 


r  -  1 8 

pAr)  =  J  —PAfoPv 


r 


(42) 


By  inserting  expressions  for  p ^  and  /?,,  as  given  by  Eqs.  (40)  and  (41).  respectively,  and  by  simplify 
ing  results  in 
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PA') 


_ K\ _ 

(. N  -  2)UK  -  N  +  1)! 


_.V  -2 


(I  +  r) 


K  *  | 


r  >  0. 


(43) 


If  we  set  q  =  1/(1  +/•),  it  is  straightforward  to  show  that 


/>,,(</)  = 


K\ 


(N  -  2)!  (K  -  N  +  1)! 


(1  ~qf 


,  K -N  +  I 


,()<</<  1  • 


(44) 


which  is  the  PDF  derived  by  Reed  et  al.,  |9|  for  the  normalized  instantaneous  S/N  p  wer  ratio  that 
results  if  the  sampled  matrix  inversion  (SMI)  algorithm  is  used.  By  substituting  q  for  r  in  Eq.  (38), 
we  obtain 


Now 


p{D  |  =  exp  - 


1 

q  a;  +  1 


riT. 


(43) 


P(D\q)  =  | 


1 


o  (K  -  N)\ 


v  exp  - 


1  + 


q  a;  +  1 


dr, 


1  + 


q  a  ‘  4-  1 


-(At  -  +  I) 


(46) 


or 


^  =  I 


K ! 


o  {N  -  2)\(K  -  N  +  1)! 


+ 


r  _s 

_s_ 

N 


+  1 


-(At  -N  +  1) 


(l-qf-2qK~N  +  ldq.  (47) 


opt 


where  a;  =  (S/N )opl .  We  note  in  Eq.  (47)  that  we  have  replaced  a;  with  (S/N)„pt,  where  (S/N)opl 
is  the  optimal  (S/N)  output  power  ratio  of  the  matched  filter  (K  =  oo).  We  can  write  this  in  this  way 
because  the  output  noise  power  of  the  Mh  channel  has  been  normalized  to  one  and  the  output  of  the 
Mh  channel  is  the  optimal  matched  filter  output. 


2.4  Probability  of  False  Alarm 

The  probability  of  false  alarm  PF  is  easily  found  from  Eq.  (47)  by  setting  (S/N)opt  =  0.  It  is 
found  that 


Pf  = 


1 _ 

(T  +  i)*  'v  +  1  ' 


An  equivalent  expression  for  PF  as  given  by  Eq.  (48)  was  also  derived  by  Refs.  I  and  4. 


(48) 
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3.0  MEAN  LEVEL  ADAPTIVE  DETECTOR  (MLAD) 

3.1  Detector  Form 

An  intuitive  form  of  adaptive  detection  is  found  by  implementing  MLAD.  Figure  1  shows  that 
the  MLAD  is  essentially  an  AMF  followed  by  an  MLD.  The  MLAD  is  designed  to  perform  detec¬ 
tions  over  a  block  of  data  by  using  just  this  block  of  data  in  determining  the  AMF  weights  and  the 
MLD  threshold.  The  MLAD  works  as  follows.  Let  there  be  A  channels  and  K  -F  1  samples  per 
channel.  Define 

x  as  the  primary  A -length  data  vector. 


xk  as  the  secondary  A -length  data  vectors,  k  =  1.2 . A\ 

Xaui!  as  the  (x|X)  -  augmented  A  X  (K  +  I)  matrix  of  input  data,  and 
*o  as  the  Xaug  X"Uj,. 


SENSORS 


r 


ADAPTIVE 

matched 

FILTER 

("=r;>> 


Fig.  I  Mean  level  adaptive  detector 
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The  A-length  weighting  vector  w  for  the  AMF  is  found  by  the  SMI  algorithm  and  is  given  by 


-  - 1 

w  =  /?0  s. 


This  weight  is  used  in  the  detection  rule  given  by 

-  H  >  K  „  H 

|w  x|2  To  £  |w  xk  t2,  (50) 

k  =  I 

Ho 


where  To  is  chosen  to  control  the  false  alarm  probability.  We  see  that  Eq.  (49)  is  the  algorithmic 
representation  of  the  AMF  and  Eq.  (50),  the  MLD. 

Note  that  we  have  included  the  primary  data  vector  in  the  R0  estimate  and,  hence,  in  the  w  esti¬ 
mate.  In  a  practical  situation,  this  might  be  done  since  it  is  more  numerically  efficient  to  compute 
one  weighting  vector  over  the  entire  data  block  than  it  is  to  compute  a  distinct  weighting  vector  for 
each  point  in  the  block.  However,  the  presence  of  the  desired  signal  (under  //,)  will  affect  detection. 
In  Eq.  (50),  the  primary  data  vector  is  varied  across  the  K  +  1  data  snapshots,  where  the  xk  used  on 
the  right  side  of  Eq.  (50)  does  not  include  a  selected  primary  data  vector. 

It  is  straightforward  to  show  that 

To  £  I  wH  x*  | 2  +  To  I  wH  x  | 2  =  T0  sH  R0  '  s.  (51) 

k  =  1 


Thus  Eq.  (50)  is  equivalent  to 

|sH  Ro'  x  [ 2  >  To 
sH  R0  1  s  <  T0  +  1 

Ho 


(52) 


where  we  note  that  0  <  7,  <  1 . 

Now  we  can  write 

R()  =  Ry  +  xxH,  (53) 


where  Rx  =  XXH.  By  using  the  matrix  inversion  lemma  [12]  and  after  some  simplification,  we  can 
show  that 
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"  -  I 


s"  /?,  '  \  |- 
1  +  \H  R ,  1  x 


1  i 


sM  R0  s  =  sH  Rx  s 


sH  /?o  X 


1  -  i 


,vH  R,  1  X 


1  +  XH  /?,  '  X 


1  -  1 


and 


xH  R0  x 


~  - 1 


xH  k; 1  x 

1  4-  xH  /?,  '  x 


(54) 


(55) 


(56) 


By  substituting  Eqs.  (54)  through  (56)  into  Eq.  (52).  the  equivalent  deteetion  rule  becomes 


> 


s "  R<  x|-  <  7*, 

Ha 


(sH  R~l  s)(l  +  xH  /?,  '  x)2  +  |  sH  Rx  ‘  x  | 2  ( 1  +  xH  R x  '  x) 


(57) 


3.2  Statistically  Equivalent  MEAD 

We  make  the  same  assumptions  on  the  input  distributions  and  statistics  as  given  for  the  GLRT 
discussed  in  section  2.  Hence,  we  can  use  the  matrix  transformations  here  that  were  used  to 
transform  the  GLRT  analysis  into  a  simpler  problem.  Specifically,  we  can  show  that  the  test  given 
by  Eq.  (57)  is  statistically  equivalent  to 


Io  (UUH)-'  v 


> 

< 

Ha 


l,1,'  (UUHV  1  1„  (1  +  vH  (UUHy[v) 


+  |  Ij  (UUH)  v  | 2  ( 1  +  vH  (UUH)  1  v)2 


(58) 


where  in  Eq.  (57)  we  have  replaced  sH  R,  '  x  with  1  o  (U(JH)  'v.  sH  Rx  '  s  with  1  ,T  (UUH)1  1(), 
and  xH  R,  x  with  xH  (UUH)  1  v.  The  quantities  U  and  v  are  defined  in  section  2.  Equations  (25) 
through  (27)  give  expressions  for  I(),  j  lj?  (UUH)  1  v|2,  and  vH  ((/f/H)  1  v.  respec¬ 

tively.  If  these  are  substituted  into  Eq.  (58),  after  somewhat  extensive  simplification,  the  following 
detection  rule  results 


H\ 


i  -  r, 

(  ■>  "3 

> 

r  2 

,  .  ZL 

< 

l  «Ti  J 

[  "ii  , 

(59) 
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We  note  for  7',  >  0.  that  if 


then  //„  is  declared. 


1  -  T, 


1  + 


_n_ 

»Ti 


<  0. 


(60) 


3.3  Probability  of  Detection 


Again  under  the  77,  hypothesis,  we  assume  that  the  desired  signal  is  a  complex,  circular  Gaus¬ 
sian  RV  (the  amplitude  is  Rayleigh  distributed).  As  in  our  analysis  of  the  GLRT. 
a;  =  a,  sH  M  " 1  s.  Assume  Eq.  (60)  is  not  true.  We  can  write  the  decision  rule  given  by  Eq.  (59) 
as 


where 


and 


«  1 2  <  T  { 


H  o 


U  ||  V  2  —  V’  j  l(  2j 


a  = 


( u] |  (a~  -F  1 )  +  ) 


2  .1/2 


T\  =  - - - ; -  7,7',. 

(q  -  T |)  (q  a-  +  1) 


(61) 


(62) 


(63) 


rj  is  defined  by  Eq.  (34),  and  q  =  (1  +  r)  '.  As  before  under  H , ,  v,  is  a  complex,  circular  Gaus¬ 
sian  RV  with  power  equal  to  'tr.  +  1.  and  u2\  is  the  same  with  power  equal  to  one.  Furthermore,  a 
is  the  same  with  power  equal  to  one.  Thus 


f0  if  q  <  T | 

P(D  |  <7,  tj)  =  5  _r 

e  otherwise. 


(64) 


The  PDFs  of  t]  and  q  are  given  by  Eqs.  (39)  and  (44).  respectively.  After  some  simplification,  we 
find  that 


p»  =  f/  i()  P{D  I 


'll 


K 1 


(N  -2)\(K  -N  +  1 )! 


I  + 


(</-7'|) 


where 


r 

r 

~\ 

s 

+  1 

q 

N 

J 

opt  j 

-  ( K  -  N  *  1 1 


(1  -q)N~2qK-s"dq,  (65) 


_S 

N 


opt 
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4.0  RESULTS 

Here,  we  present  results  on  the  detection  probability  Pn  for  the  GLRT  and  MLAD  \s  the 
independent  parameters:  the  probability  of  false  alarm  P,  ;  the  steady  state  (K  =  oo>  S  N  output 
power  ratio  ot  the  matched  filter  (S/N)opt:  the  number  of  independent  samples  per  channel  K  of  the 
secondary  data;  and  the  number  of  input  channels  N.  We  set  K  =  MN,  where  M  is  a  positive 
integer,  and  M  is  an  independent  parameter  called  the  degrees  of  freedom  (DOF)  factor.  Figures  2 
through  17  show  plots  of  Pn  vs  (S/N)opt  for  PF  =  10  6.  10  m.  N  =  2.  5.  10.  30.  and  the  two 
detection  schemes. 


S/N  OPT  (dB) 


Fig.  2  -  GLRT:  Pn  for  N  =  2.  Ph  =  l  D-fi 
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S/N  OPT  (dB) 


Fig.  5  -  GLRT:  P„  for  N  =  30.  P h  =  I  D-6 


S/N  OPT  (dB) 


Fig  6  -  G',RT:  Pn  for  ,V  =  2.  Ph  =  I  D- 10 
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S/N  OPT  (dB) 


Fig  7  -  Cil.RT:  P„  lor  ,V  =  5.  Ph  =  !  D-IO 


S/N  OPT  (dB) 


Fig.  8  -  GLRT :  Pn  for  N  =  10,  Pv  =  I.D-IO 
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Fig.  9  -  GLRT:  Pn  for  N  =  30.  Ph  =  l  .D-10 
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Fig  10  -  MLAD:  Pn  for  N  =  2.  P,  =  l.D-6 
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Fig.  13  -  MLAD:  P„  for  N  =  30.  Ph  =  l.D-6 


S/N  OPT  (dB) 


Fig.  14  -  MLAD:  P„  for  N  =  2.  P,  =  I  D- 10 
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S/N  OPT  (dB) 


Fig.  15  -  MLAD:  Pn  for  N  =  5.  Ph  =  l.D-10 
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Fig.  16  -  MLAD:  P„  for  ,V  =  10.  =  l  .D-10 
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Fig.  17  -  MLAD:  P„  for  N  =  30,  Py  =  I  D- 10 


Kelly  [1|  defines  the  (S/N)  loss  of  an  adaptive  detector  as  the  difference  required  to  obtain  a 
given  Pn  between  a  steady-state  (M  =  oo)  detector  and  the  transient  state  with  all  other  independent 
parameters  being  equal.  We  note  that  as  K  —  oo  (or  M  —  oo),  the  GLRT  and  MLAD  are  identical 
and,  hence,  their  PD  performance  is  identical.  We  define  M3dB  to  be  the  DOF  factor  such  that  the 
(S/N)  loss  is  nearest  to  3  dB.  For  both  the  GLRT  and  MLAD,  we  make  the  following  observations 
from  Figs.  2  through  17. 

(1)  Both  detectors  are  slower  to  converge  to  their  optimal  value  ( M  =  oo)  for  the  lower 

ordered  matched  filters.  For  example,  for  the  GLRT  for  most 

PDs  (0.1  -  0.9).  PF  =  10~6,  if  N  =  2,  then  A/3dB  =  6,  if  N  =30.  M3dB  =  2. 

(2)  There  is  a  diminishing  return  in  convergence  by  using  a  larger  DOF  factor,  i.e.,  there  is 
more  improvement  in  performance  in  going  from  M  to  M  +  1  than  M  +  1  to  M  +  2. 

(3)  Convergence  slows  for  decreasing  PF.  For  example,  for  the  MLAD  for  most 

P[)  (0,1  -  0.9)  and  N  =  5,  if  PF  =  10~6,  then  M3dB  =  5;  if  PF  =  10' 10 ,  then 

~  6- 

We  note  that  these  trends  were  also  observed  by  Kelly  [1J  for  the  GLRT,  where  the  signal  is  modeled 
as  having  a  constant  amplitude. 
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In  comparing  the  GLRT  with  the  MLAD,  we  make  the  following  observations  from  Figs.  2 
through  17: 

(4)  The  P[,  performance  for  fhe  GLRT  is  always  superior  to  the  PD  performance  of  the 
MLAD  with  all  independent  parameters  being  equal  between  the  two  detectors. 

(5)  For  low  DOF  factor  (M  =  2,  3).  the  GLRT  has  much  better  PD  performance  than  the 
MLAD  with  all  other  independent  parameters  being  equal  between  the  two  detectors.  For 
example,  if  N  =  2,  Pt  =  10“6,  (S/N)opt  =  30  dB,  and  M  -  2  for  the  two  detectors, 
then  for  the  GLRT,  PD  =  0.7  and  for  the  MLAD,  PD  =  0.09. 

(6)  The  A/,dB  fcr  the  GLRT  and  MLAD  are  approximately  the  same  with  the  greatest  differ¬ 
ence  occurring  for  low  DOF  factors  (N  =  2).  For  example,  for  most 
PDs  (0.1  -  0.9),  PF  =  10~6.  if  N  =  2,  then  MjdB  =  6  for  the  GLRT  and  B  =  7 
for  the  MLAD.  For  the  same  case  but  N  =  30.  the  A/,dB  is  identical  for  both  the  GLRT 
and  MLAD. 

^  H 

The  GLRT  outperforms  the  MLAD  because  the  MLAD's  MLD  samples  w  x*.  k  =  1,  2,  ...  K  are 
being  contaminated  by  the  desired  signal  in  x  since  w  is  computed  by  using  a  sampled  covariance 
matrix  RQ,  which  contains  the  desired  signal  in  some  form.  This  results  in  an  increased  average 
power  in  the  MLD’s  samples,  which  results  in  an  increase  i"  f>ie  MLD’s  threshold,  which  results  in 
fewer  detections.  The  increase  in  the  residue  power  in  ♦he  rang  :  cells  (for  the  radar  problem)  around 
the  primary  range  cell  because  of  signal  co"G  .mation  was  also  reported  in  Lewis  and  Kretschmer 
[131  and  Ger/ach  (I4j. 

Of  the  three  comparisons  between  the  CLAY  and  .he  MLAD,  observation  (6)  above  is  possibly 
the  most  significant.  As  we  will  show  in  section  5,  the  GLRT  requires  more  numerical  computions  to 
implement  than  the  MLAD.  However,  as  we  note  from  the  observation  (6).  the  required  number  of 
samples  will  be  approximately  the  same  for  both  detector  schemes  in  order  to  have  only  a  3  dB  S/N 
detection  loss. 

5.0  COMPUTATIONAL  CONSIDERATIONS 

Here  we  briefly  compare  the  computational  requirements  for  the  MLAD  and  GLRT  detection 
schemes.  First  we  consider  the  GLRT.  We  modify  the  form  given  by  Eq.  (1).  Now 

Rx  '  =  (i?0  -  x  xV 


-  - 1  #0  1  X  XH  Rn  1 

*0  + 


'0 


1  -  xH  R 


i 


(68) 


o 


If  the  expression  on  the  right-hand  side  of  Eq.  (68)  is  substituted  into  Eq.  (I)  for  Rx  ',  after  some 
simplification,  the  GLRT  detection  rule  becomes 


u  "  ,,  11  I 

SH  *0  X  1  ~  _ 1 _  > 

SH  Rn  '  S  (1  -  xH  R()  1  x)  u 


(69) 
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where  t'  =  //( 1  -/).  It  we  compare  this  equivalent  form  of  the  GLRT  to  the  MLAD  form  as  given 
by  Eq.  (52).  we  see  that  the  GLRT  test  statistic  is  merely  the  MLAD  test  statistic  multiplied  by 

(1  -  \H  R0  x)-1.  It  is  straightforward  to  show  that  0  <  1  -  xH  /?„  x  <  1,  and  thus  the  GLRT's 
test  statistic  is  always  greater  than  the  MLAD’s  test  statistic  under  H0  or  //, . 

We  see  that  it  is  the  additional  computation  of  1  -  xH  /?„  x  for  the  GLRT  that  is  the  major 
computational  difference  between  the  two  detection  schemes.  There  are  many  schemes  for  imple¬ 
menting  the  adaptive  matched  filter  in  open-loop  fashion:  SMI,  Gram-Schmidt,  Givens  rotation. 
Householder,  etc.  combined  with  a  block  processing  or  sliding-window  technique  (of  which  there  are 
many  variations).  No  efficient  method  has  been  proposed  in  implementing  the  GLRT,  although  a 
variant  that  uses  the  normalized,  fast  orthogonalization  network  described  in  Refs.  15  and  16  could  be 

used  to  compute  R{)  x  and,  thus,  xH  R0  x.  However  the  implementation  complexity  of  this 
scheme  is  significantly  greater  than  that  of  the  AMF.  The  efficient  implementation  of  the  GLRT  is 
left  as  a  topic  of  future  research. 

6.0  SUMMARY 

Two  schemes  for  adaptive  detection  have  been  analyzed  and  compared:  Kelly's  GLRT  and  the 
MLAD.  Detection  performance  Pn  was  predicted  for  the  two  schemes  under  the  assumptions  that  the 
input  noises  were  Gaussian  random  variables  that  were  temporally  independent  but  spatially  corre¬ 
lated,  and  the  desired  signal  was  Rayleigh  distributed.  PD  was  computed  as  a  function  of  the  false 
alarm  probability,  the  number  of  input  channels,  the  number  of  independent  samples-per-channel,  and 
the  matched  filtered  output  S/N  power  ratio. 

In  this  analysis,  the  GLRT  was  shown  to  have  better  detection  performance  than  the  MLAD. 
The  difference  in  detection  performance  increases  as  one  uses  fewer  input  samples.  However,  the 
require^  -.lumber  of  samples  necessary  to  have  only  a  3  dB  detection  loss  for  both  detection  schemes 
is  approximately  the  same.  This  is  significant  since  for  the  present,  the  MLAD  is  considerably  less 
complex  to  implement  than  the  GLRT. 
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